ESCAPE RATES FORMULAE AND METASTABILITY FOR 
RANDOMLY PERTURBED MAPS 



WAEL BAHSOUN AND SANDRO VAIENTI 

Abstract. We provide escape rates formulae for piecewise expanding interval 
maps with 'random holes'. Then we obtain rigorous approximations of invari- 
ant densities of randomly perturbed metabstable interval maps. We show 
that our escape rates formulae can be used to approximate limits of invariant 
densities of randomly perturbed metastable systems. 



1. Introduction 

A dynamical system is called open if there is a subset in the phase space, called a 
hole, such that whenever an orbit lands in it, the dynamics of this obit is terminated. 
In open dynamical systems, long-term statistics are described by a conditionally 
invariant measure and its related escape rate, measuring the mass lost from the 
system per unit time [T31 [H] . For the past three years there has been a consider- 
able interest in describing the escape rate of an open system as a function of the 
hole's position and size. In [8j it was observed that for the doubling map, given two 
holes, Hi and H2 of the same fixed size, each centred around a periodic point, say 
X\,X2 respectively, with the period of x\ smaller than that of x 2 , then the escape 
rate through H\ is bigger than that through H2. Later in [57] it was shown that, 
by shrinking the hole to a point, say xo, the escape rate depends on two things: i) 
whether x is periodic or not; ii) the invariant density of the corresponding closed 
system. Following the success of [H EJJ , other researchers studied this phenomenon 
for different types of systems [TJ HH ITTQ These results have lead to insights in 
studying metastable dynamical systems which behave approximately like a collec- 
tion of open systems: the infrequent transitions between almost invariant regions in 
a metastable system are similar to infrequent escapes from associated open systems 

[sang iniiii]. 

This research topic is currently very active in ergodic theory and dynamical 
systems. In part, this is due to the interesting applications of open and metastable 
dynamical systems in physical sciences. For instance, open dynamical systems are 
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used to study transport in heat conduction 21 . They also play an important 
in astronomy [33] . In molecular dynamics, almost invariant regions of metastable 
systems are used to identify sets where stable molecular conformations occur [31 j . 
In astrodynamics, metastable systems are used to recognize regions from which 
asteroids escape is infrequent [11] . Moreover, metastable dynamical systems have 
been recently used to develop realistic models of atmospheric and ocean circulation 
[TU1 [THl ED] ■ In these models, it is believed that metastable states^lie behind long 
term global circulation patterns, and form large scale barriers to transport. 

In this paper we first study piecewise smooth and expanding interval maps with 
'random holes'. We prove that the random open system admits an absolutely con- 
tinuous conditionally stationary measure (accsm) . Moreover we obtain escape rates 
formulae (first order approximations) depending on the position of the holes. An 
older escape rate formula for randomly perturbed maps with specific holes positions 
and a specific distribution on the noise space was obtained in [S]. Our result is in 
the spirt of the recent trend of describing the escape rate as a function of the po- 
sition of the hole (TJ [U [T7J [57J . In particular, our result generalizes Keller-Liverani 
escape rates formulae 27 which were obtained for deterministic perturbations. 

In the second part of the paper we study random perturbations of interval maps 
that initially admit exactly two invariant ergodic densities. However, under random 
perturbations that generate random holes and allow leakage of mass between the 
two initially ergodic subsystems, the random system admits a unique stationary 
density. We show that such a density for the metastable random system can be 
approximated in the L 1 norm (with respect to the ambient Lebesgue measure m), 
by a particular convex combination of the two invariant ergodic densities of the 
initial system. In particular, we show that the ratio of the weights in the convex 
combination is equal to the ratio of the averages of the measures of the left and 
right random holes. Moreover, as a by product of our escape rates formulae, we 
show that these weights can be also identified as the ratio of the escape rates from 
the left and right random open systems. 

We finish the introduction by addressing possible generalizations of our work. We 
foresee two natural generalizations. The first generalization concerns the extension 
to piecewise expanding maps in higher dimensions. Two difficulties arise in this 
setting. First, obtaining a Lasota-Yorke inequality for the operator associated with 
the open system with random holes, see Lemma |2.6| would require a more careful 
control of the variation, whatever the Banach space B that operator acts on is, 
of the characteristic function along the boundaries of the partitions, which are 
now codimension-1 piecewise smooth submanifolds and not points as in the setting 
of this paper. The second difficulty is to get local smoothness of elements of B 
around the holes similar to Lemma 1 of [55] . We believe that the latter property is 
fundamental to obtain results on approximating invariant densities of metastable 
systems, even in the one-dimensional case. 

The second generalization is the extension to non-uniformly expanding maps. 
Mctastability for deterministic maps with indifferent fixed points was treated by 



Metastable states correspond to eigenvalues of a transfer operator which are very close to 1, 
a phenomenon which is very transparent in the metastable systems of this paper. 
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us in |3] ■ One of the main tools we used there was to reduce the study of metasta- 
bility to an induced subsystem where the map is uniformly expanding, and then 
pullback the result to the original system. It is not clear how this technique can be 
implemented in the presence of noise. A few recent results [2J [35] used induction 
to achieve stochastic stability in some classes of one-dimensional maps. It would 
be interesting to explore the possibility of apply those techniques to metastability. 
Supposing the analogous of our Theorem |3.4| can be proved in the random case, 



it would be also interesting to explore if something similar to Corollary 3.5 (see 
section 3.3) would still hold. In this paper, the latter is related to the exponen- 
tial escape rate which is in turn associated to the existence of a spectral gap of 
the transfer operator. This is a 'luxury' that is not available in a non-uniformly 
hyperbolic setting. 

In section [2] we obtain escape rates formulae for expanding interval maps with 
random holes. Our main result in this section is Theorem 12.11 In section 13.31 we 
study rigorous approximations in the L 1 -norm of invariant densities of randomly 
perturbed metastable maps. Our main result in this section is Theorem |3.4| 



2. Escape Rate Formulae For Randomly Perturbed Maps 

2.1. Notation. Throughout the paper we use the following notation: (I,Q3,m) is 
the measure space where I — [0,1], 93 is the Borel c-algebra and m is Lebesgue 
measure. For / € 53,771), we define 

Vf = inf{var/ : / = / a.e.}, 

where 

t-i 

varf = sup{J^ \f(x j+ i) - f(xj)\ : = x < x x < ■ ■ ■ < x t = 1}. 

3=0 

We denote by BV{I) the space of functions of bounded variation on / equipped 
with the norm || • \\bv = V(«) + || • ||i, where || • ||i is the L 1 norm with respect to 
m . 

2.2. Piecewise expanding maps. Let T : I — >• / be a map which satisfies the 
following conditions: 

(01) There exists a partition of /, which consists of intervals {Ii}f =1 , U fl Ij = 
for i 7^ j, Jj := [c i , c i+1 ] and there exists S > such that T i :— T\( c . 0;Ci+1 ) is 
C 2 which extends to a C 2 function T^o on a neighbourhood [c^o ~ 5, c^+i o + S] of 

h ; 

(02) inf BeA c \T'(x)\ > k- 1 > 1, where C = Ko}? = r 

(03) T preserves a unique acinj^]/i which is equivalent to m. Moreover, the system 
(I, /i, T) is mixing. 

2.3. Random Holes. Let (wfc)fegN be an i.i.d. stochastic process with values in 
the interval fl £ = [0, e], e > 0, and with probability distribution 9 £ . We fix z g (0, 1) 
and we associate with each to £ f2 e an interval H u such that z € H u C H e , and we 
assume that H e C H £ > for e < e' . Further, 

(04) We assume that T is continuous at z; this assumption will be explicitly used 
in the proof of Theorem |2.1| In order to apply the results of this section to section 



The existence of fi follows from the well known result of 1281 . 
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|3.3[ we will also assume that the density p is continuous at z. 




Figure 1. A map with random holes H ull and H U2 sitting inside 
H e . Note that H Ul and H^ 2 are not nested. 



2.4. Transfer operator of the random open system. Our goal in this section 
is to study the existence of absolutely continuous conditionally stationary measures 
(accsm) and their associated escape rates through the random holes H M around 
a given point z under the dynamics of T. An example of an interval map with 
random holes is shown in Figure 1. We set 

:— I \ H^, 

and define for / G *8, m) 

(2.1) C e f(x) := / C (fl x J(x)d9 s , 

where Co is the transfer operator (Perron-Frobenius) (4j [Tj associated with T; i.e., 
for / G L 1 (7,03, to) we have 

Cof{x) = lk M ww\- 

y:=T 1 (x) 

The transfer operator C s will be used to prove that the random open system admits 
an accsm with exponential escape rate. By using the fact that for any measurable 
set A and integrable function / 

Ia^o/ = Co{l T -x A f), 

we obtain 
(2.2) 

£"f{x) = [ ■ ■ I ^(/Ix,,, nr-^n-nr-c-DA,. )(i)i»s(«i) ■ • • * e (wn) 

def 



£o(/ 1 x J1 nT- 1 x^n---nT-(™- 1 )x^„ )( x ) d ^T ( w ) 
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where Of :— ngL^e and u := (u>i, uj 2 , . ..)• 

2.5. Statement of the main result of section [2j The following theorem pro- 
vides random versions of Keller-Liverani escape rates formulae [27] . 

Theorem 2.1. For sufficiently small e > 0, there exists an e e , < e e < 1, and a 
g G BV(I), g £ (x) > 0, with f I g £ .dm = 1, such that 

^e9e = e £ g £ . 

Moreover, 

(1) If z is a non-periodic point of T , then 

lim -= ; — ; — - = 1. 

(2) If z is a periodic point of T of minimal period p, T v is C 1 in a neighbour- 
hood of z, and the probability distribution 6 £ on the noise space Q £ satisfies 
the following condition: 3 v > 1 such that 

(C) c {uj- m(H u ) G (e - e v , e)} > 1 - e v ; 
then 

1 - e e 1 l 

lim -= ; ; .„ , r = l 



>o v{Hj)M e (u>) \(TP)'(z)\- 
Remark 2.2. 



(1) In (2) of Theorem 2.1 condition (C) means that for a fixed s > most of 
the random holes are not too small when compared to the size of H e . 
Notice that the two "u" in the condition can be different, provided they 
remain larger than l. We considered the same v for simplicity Condition 
(C) is needed to insure a lower bound on the escape rate in case (2) of 
Theorem O 

(2) For the existence of the accsm we required e to be sufficiently small since 
we use the perturbation result of |26j in the proof our Lemma 2.7 The 
strict positivity of the density g £ also follows from [26] since /i is equivalent 
to m. 

Remark 2.3. Since C £ g e = e £ g e , with Jj g £ dm = l, by using the definition of C £ , it 
follows that: 



/ d9 £ (uj) / g e l x ^,dm. 
Jn e Ji 



Set v £ to be the Borel probability measure: 



v e (A) := — I I l A ~Lx^g £ dQ £ dm, 



-e J I JO, 



It is very easy to construct examples that satisfy condition (C). For instance, suppose we 
associate to any ui £ H e a symmetric hole around the point z. Then: (i) H e = [z — §; z + |]; 

= [z — ^; z + tt]; (ii) m(H s ) = e; m(/f U j) = u. Let us choose an absolutely continuous 9 e 
with density d £ . Our condition (C) will be satisfied whenever f^_ e v d £ (ui)dw > 1 — e v which can 
be obtained, for instance, by taking d s = e~ v on the interval (e — e v ,e) and everywhere else. 
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where A C I is a measurable set. By using the fact that C e g £ = e £ g e and the 
definition of £ e , we obtain that v € satisfies the following: 

(2.3) u e (A) = — f ^(wji/^rUnrX); 



(2.4) e e = / ^(^^(T" 1 ^); 
and 

(2.5) e n £ v e {A) = I d6'?{uj)v £ (T- n AnT~ 1 X UJl C\T- 2 X ul2 n---C\T- n Xu n ). 

Definition 2.4. We will call — lne e the escape rate for the random system, and v t 
the absolutely continuous conditionally stationary measure. 



Remark 2.5. It is interesting to remark that the items (1) and (2) of Theorem 2.1 
provide explicit perturbation formulae for the leading eigenvalue of the operator 
C e . This operator enjoys all the properties of what Keller recently called rare 
events Perron- Frobenius operators. 



2.6. Proof of Theorem 2.1 , We start this subsection by proving a uni form 
Lasota-Yorke inequality for C e and £ . A deterministic version of Lemma 2.6 
can be found in [2"5] . 

Lemma 2.6. There exists a 7 € (0, 1) and constants A,B>0 such that for any 
n > 1 and f G _BV A (/) we have 



WCUWbv < Ar\\f\\BV + B 

\\CU\\bv < Ar\\f\\BV + B 

Proof. We begin to compute the total variation, leaving the L 1 estimate at the end. 
Let Z n = Z V T~ X Z • • • V T~( n ~^Z, where Z = {Ii}f =1 . Let J„ = p^yy ■ For 
A e 2" and / e BV(I): there is a k € (0, 1), such that 

(2.6) Vl T « A {J n f) o T-» < 2n n V A f + 1 --r f \f\dm. 

In particular, 

(2.7) V r £S(J)< J2 Vl T r>A(J n f)°T,- n <2 K n V I f+-. 1 [ \f\dm. 

Now, let us consider a fixed random path of length n, (u>i, uj%, . . . ,oj n ), f € BV{I) 
and define 

■9 := f^-X u , nr- 1 X„,n-nT-("- 1 )X„„ 



and observe that g S -BF(I). Then using (2.6 1 and noticing that 

n A 



(x Ul n t~ 1 x UJ2 n • • • n T-^-^X^ 
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consists of at most n + 1 connected components, we obtain 
VlT"-A(Jnflx ul nT- 1 X aa n-nT-l»- 1 )x Un ) ° T\a 

= Vl T ~A(J n g) o T,-" < 2 K n V A g + -. 1 ^ / |g|rfm 

|yl mm AeZn m(A) J A 

<2 K "(V A / + 2(n + l)sup|/|) + ^ 1 — / |/|dm 

s6 A mm AeZ n m{A) J A 

<2K n (2n + 3)V A f+— 1 ^(2«; n (2n + 2) + l) / \f\dm. 

miiiAeZ" m(A) J A 



Summing over A G Z" in (2.8 1 we obtain 
(2.9) n 

Vl^O nT~ 1 x 1Jo n---nT-(™- 1 )x 0J „ ) 



< 2 K n (2n + 3)y//+ — v (2K™(2n + 2) + l) / |/|dm 

mm AeZ i m{A) 



Since the inequality in (2.9) does not depend on lu and 



/ •■■/ -Co(/ 1 x tt1 nT-iX W9 n-nT-(«-i)X w X^d^^) • • • d(9 £ (w„), 

we also have 

(2.10) Vr£™/ < 2 K "(2n + 3)V>/ + — - r -(2 K "(2n + 2) + 1) / \f\dm. 

mm AeZ « m(A) Jj 

The estimate on ||£"/||i is easy. Indeed, this can be done by splitting / into its 
positive and negative parts and by using the linearity of the transfer operator. 
Therefore, we may suppose that / is non negative. This allows us to interchange 
the integrals w.r.t. the Lebesgue measure and 6*^° and to use duality. In conclusion 
we get 

< J Ifllx^nT-^x^n-nT-^-nx^Jdm^ 

Since there exists no and 7 € (ft, 1) such that 2ft n °(2no + 3) < 7" , we c an cho ose 
A := 2n + 3, B := sup„<„ minA J n m(A) 2(n + 1)^ and use ( |2T0| , Q to 

obtain a uniform the Lasota-Yorke inequality for £ e and Cq. 

□ 

Lemma 2.7. Consider L E : BV(I) -t BV(I). Then there are e e £ (0,1), <^ 6 € 
£?V(7), a probability Borel measure v s and linear operators Q e : BV(I) — > BV{I) 
such that 

(1) eJ 1 £ e = (p e ®v e + Q e ; 

(2) C £ ip e = e e ip e , v £ C e = e e v £ , Q £ ip e = 0, v e Q e = 0; 

(3) E^=o su Pe \\Q1 1 \\bv(i) < 00; 

(4) m{ip £ ) = 1 and sup E ||<£ E ||w(/) < °o; 

(5) 3 C > sucft tfiai 

r^e := sup I / (£0 — £ E )ipdm\ — !• as e — >• 0, 
II</>IIbv(j)<1 Ji 

and 

r? e -||(£o-4)p||w(7) <C|A £ |, 
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where A £ :— m((Lo — £-e)p)- 
Proof. We have for any -0 G BV(I): 

(2.11) (C - C e )ip 

Since tp is also in L°°(m) we have 
(2.12) 



sup 

II^IIbv(d<i •" 



(Co - £ e )V)rfm| < / m(H u )d0 £ < m(H e ) ->• as e -)• 0. 



Thus, by Lemma 2.6 and (2.12), using the abstract perturbation result on the 
stability of spectrum of transfer operators [2E], we obtain (l)-(4) of the lemma. 
Now by (2.11 ), notice that 



(2.13) 

and 
(2.14) 



% < A e := [ m(H u )d6 £ , 
Jsi, 



p(H u )d0 £ 



Moreover, using (2.6), condition (03), and calling: (i) K = max(A7,i?); (ii) 
G = Z\\p\\bv\ (hi) p m := infsei p(x),x — m a.e., we obtain 

||(£ -£ s )p\\bv(i) = II / Co(pl H J)de e \\ B v{i) 



(2.15) 



<K [ (V(pl H J + \\pl H J\i)d9 £ < 



KC / n « ^( H ^)d0e 

Pm. A F 



Thus, by (2.13), (2.14) and (2.15), 



Ve ■ 11(A) ~ C £ )p\\bv(i) < 



KC 



□ 



We are now ready to apply the abstract perturbation result of 27 j. Indeed, 
Lemma 2.7 shows that all the conditions imposed in [27 are satisfied by our random 
systems. 

Lemma 2.8. 

(1) For sufficiently small e, A £ ^ 0. 

(2) If for each integer k > the following limit: 

m((Co-£e)C k e (L -£ e )(p)) 



q k := lim q k , £ := lim 

e->-0 e^O 



A, 



exists, then 



lim — t — - 



fe=0 



Proof. By condition (03), for sufficiently small e > 0, A e > 0. Thus, Part (2) of 
the lemma follows from the abstract perturbation result of [S7J , since we will show 
in the proof of Theorem 2.1 that for each k > the limit in (2) exists. □ 
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Proof. (Proof of Theorem 2.1 ) We start by computing qk, £ and show that the limit 
in (2) of Lemma 2.8 holds. Throughout the computation in (2.16) below, we will 
use repeatedly the fact that for any measurable set A and integrable function / we 
have IaCqI = C (l T -i A f ). 

Consequently, by (2.111 and (2.2 1, we obtain 

(2.16) 

m((£ -4)4 fe (A>- A)(p)) 



£l E J I 



n E J i 



n E J i 



Q e J I JO. 



((£„ - C £ )C k JC - C £ ){p))dmd9 £ {u) 

£ (1 H „ ' £e(£o - £ s ){p))dmde e (u) 

1 H „ ■ C k £ (C Q - Ce)(p))dmd6 e (uj) 

1 h w ■ -Cq1x w ,n-nT-c*-i)x M ■ {£-o - t £ )(p))d6 e {u))dmd6 e {w) 

Cq(1 t -i, Hw ■ l Xuin ... nT -(k-i)x Uk ■ (A - C £ ){p))d9 £ (Q)dmd8 £ (uj) 

ir-KH^ ■ ii^n-nr-c-nx^ ' (A) - C e )(p))dmd9 e (u))d6 e (u>) 

1t-*h„, ■ ix^n-nr-c-nx^ ■ (C lH u ,p)dmdO e (Q)d6 e (u)d9 e (u') 

/ / / / £oQ-t-v>+vh u ■ 1 T-ix„,n-nT-i'X u • lH u ,p)dmd9 e (u)d6 e (u})d6 e (u}') 
Jn E Jn E J i Jn E 

n {[t-^h^ n T~ 1 x Ul n • • • n t^x^} n d6 e {w')d6 e (u))d6 e {uj). 



n E J i Jn 



n E Jn E J i 



n E Jn E J i Jo. 



n, Jn E Jn 



Now, if z is not a periodic point of T, since T is continuous at z, for sufficiently 
small e > 0, 

{T~ 1 x UJ1 n • ■ ■ n T~ k x Uh n t~^h u ] n h u , = 0. 



Therefore, (2.161 implies qu, £ = qk = 0. This proves (1) of the theorem. 



To prove (2) of the theorem, since z is a periodic point of minimal period p and 
the density p is essentially bounded from below, the quantity we have to compute 
reduces to 

JJm(T-PH u nH u ,)d0 e (u;)(W e (u') 
{ ' f m{H„)dB € {w) 

Since T p is continuous around z we can suppose that, by taking e small enough, all 
the H^, \uj\ < e, will be contained in a neighborhood of z where T p is continuous 
and monotone. This in particular implies the following fact that we will use later 
on: if we call H UtP the unique connected component of T~ P H UJ which contains z, 
it will be the only connected component of T~ P H U1 which intersects any other H u i. 
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By the differentiability assumption on T p and the local change of variable we have 
m(T- p H u nH u ,)= f dm= f \DT p (y)r 1 dm(y) 

= [ iDTPiy^dmiy). 

JH^nTPH,,, 



Equation (2.18) gives immediately the upper bound 



(2.19) m{T-' p H UJ nH u ,)< sup \DTP\- 1 m{H u ) < sup \DTP\~ 1 m(H u ). 



By (2.191 and the continuity of DT P at z we obtain 



e^o J m(H u )d6 € {ijj) 

For the lower bound, since T p is expanding, we have 

m(T"^n^)> / \DT p {y)\- 1 dm{y) 
(2.21) =/ \nF>{y)\- l dm{y) - [ \DT p (y)^ 1 dm{y) 

> inf \DTP\- 1 m{H u ) - [ \DT p (y)\~ 1 dm(y). 



Therefore, 
(2.22) 

J f mjT-PH^ n H u ,)d8 e (Lj)d6 e {uj') 

J m{H u )d6 t {bj) ~ J J H* 

1 1 Ih„/h„, \ DTP (y)l^ 1 dm(y)dd € (uj)d8 e (uj') 
J m(H u )dB e (u) 



> / / inf | DT p \- l d0t(w)d0t(u') 



Let us consider the second expression on the right hand side of (2.22 1 and show 
that it vanishes when e goes to zero. First, we have 

/ \DT p (y)\~ 1 dm(y) < sup \DT* I" 1 m(H u AH u ,). 

Now, let us call G w the complement of H u in H £ ; we immediately have by the 
condition (C): m(H u n H u >) = m(H u >) - m(H u , n G u ) > (e - e v ) - e v > e - 
2e v . This implies that m(H u /H u >) = m(H u ) - m(H u n H u >) < 2e v and therefore 
m(H u AH u ,) < Ae v . Let us call F £ := {uj G Vt e - m{H^) G (s - e v ,e)}; Q £ := 
{{oj,u') eS! E x rt e ; m(H u AH u >)} < 4e v . Notice that Q e D F e x F £ and recall that 
O e {F B )>\-e v . 
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Then 



Sm(H u )d6 e {u) 



< 



< 



1 



< 



< 



/ m(H u )d6 e (u;) j jQc 

m(H u AH u ,)d6 e (u)de e (u>')} 

' [4e* + e 2 (F £ c ))] 



m(H u AH ul )d0 e (u)d0 e (cj') 



1 



[4e« + < 



4e u + e 



J^miH^dOM J (e-e-)(l-e-) 

which goes to zero when e tends to 0. Consequently, taking the limit as e goes to 
in (2.22), by the continuity of DT P at z, we obtain 

JJ miT-PH^ n H u ,)d6 e (u)dd e (cj') 



(2.23) 



lim 

e->0 



/ m{H u )de e {u>) 



> |£)T P | _1 (2) 



Thus (2) of the theorem follows by (2.20) and (2.23). 



□ 



3. Metastability of randomly perturbed maps 

In this section we study random perturbations of Lasota-Yorke maps. We assume 
that the initial system admits exactly two ergodic invariant densities. Then under 
random perturbations which allow leakage of mass through random holes, we will 
show that the system admits a unique stationary density. Our goal is to show 
that such a density can be approximated in the Z^-norm by a particular convex 
combination of the two invariant densities of the initial system. A deterministic 
Lasota-Yorke system of this type was studied in [35]. Deterministic intermittent 
systems of this type were studied in [3] . Other recent results on metastable random 
dynamical systems have been obtained in [20l [24] . We will establish a close link 
between the escape rate formulae which we obtained in the previous section and 
invariant densities of randomly perturbed metastable systems. We first introduce 
the class of maps of this section. 

3.1. The initial system. Let T : I — s- I be a map which satisfies the following 
conditions: 

(Al) There exists a partition of I, which consists of intervals {Ii}J =1 , 1% fl Ij = 
for i ^ j, Ii :— [ci t o, Q+i.o] and there exists S > such that T^o := T\/ Ci 0;Ci+1 ) is 
C 2 which extends to a C 2 function T^.o on a neighbourhood [q.o — 6, + S] of 

h ; 

(A2) inf KeI \ Co \T'(x)\ > > 2, where C = {o,o}Ui- 

(A3) 3 6 in the interior of I such that T\i t C /«. where * G {I, r}, is an interval 
such that Ii U I r = I and Ii D I r — {b}. 

(A4) Let H := T _1 {&} \ {b}. We call H the set of infinitesimal holes and we 
assume that for every n > 1, (T"C ) D H Q = 0. 
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(A5) T admits exactly two ergodic a.c.i.ms fj,*, such that each /i* is supported on 
I* and the corresponding density p + is positive at each of the points of H a n 



Remark 3.1. As shown in [22], assumption (A4) guarantees that the invariant 
densities of the two ergodic measures are continuous at each of the infinitesimal 
holes. 

Assumption (A2), more precisely the fact that Pq 1 is strictly bigger than 2 instead 
of 1, is sufficient to get the uniform Lasota-Yorke inequality of Proposition 3.2 as 
is explained in section 4.2 of [22]. Finally, we consider T to be bi- valued at the 
points Cj o of discontinuity in Cq and we take T(c.j i o) be both values obtained as x 
tends to c^o from either side, and T(cj ; o±) will be the corresponding right and left 
limits. 



An example of a map To satisfying the above assumption is shown in Figure 2. 
We denote the transfer operator (Perron- Frobenius) [U[7] associated with T by Cq. 




Figure 2. A typical example of the initial system T. 



3.2. Random perturbations on a finite noise space. As in the previous sec- 
tion, let (wfe)fcgN be an i.i.d. stochastic process with values in £l £ and with proba- 
bility distribution 9 e . However, we assume in this section that £l e = {so, Si, . . . , s^,} 
is a finite set, and 9 e is an atomic measure on ilj^j We associate with each we!! E 
a map T u : I — > J with T = T So , T e = T SL , and we consider the random orbit 
starting from the point x and generated by the realization cu n = (wi, u>2, ■ ■ ■ ,w n ), 



In fact, the assumption that the set fl £ is a finite set is not need in some of the proofs of 
this section. In particular in Lemmas |3.7| and |3.8| However, this assumption is needed to be 
able to identify the jumps in the stationary density p e of the random system. It seems that 
the assumption of this section that £7 e is a finite set is the best one can do when dealing with 
approximating the invariant density of the random metastable system. See Remark |3.10| where 
we explain why the problem becomes intractable if one relaxes this assumption even to the case 
where f2 e is countable but not finite. 
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defined as : :— T LUn o ■ ■ ■ o T Ul (x) (random transformations). This defines a 
Markov process T £ with transition function 

V(x,A) = f l A (T u {x))dB B (u), 

where A € 55(7), a; € J and 1a is the indicator function of a set A. The transition 
function induces an operator U* which acts on measures A on (7, 03(7)) as: 

U*A(A) = [ [ l A (T u {x))d6 e (uj)dA(x) = f U £ l A (x)dfi £ (x), 
JiJn e J i 

where 1l £ is the random evolution operator acting on 7^ functions g: 

(3.1) U e g = [ goT u d6 e (uj). 

A measure p £ on (7, 03(7)) is called a T £ -stationary measure if and only if, for any 
A e 05(7), 

(3.2) Wp £ (A) = fx e (A). 

We are interested in studying the metastability of T e -stationary measures which 
are absolutely continuous with respect to m: let us call them acsm. By ( |3.1[ ), one 
can define the transfer operator C £ (Perron-Frobenius) acting on 7 1 (7, *B(7), m) by: 

(3.3) (C £ f)(x) = [ C u f(x)d9 E (w), 
which satisfies the duality condition 

(3.4) J gCJdm = Ju e gfdm 

where g is in some subset of 7^ and £ w is the transfer operator associated with T u . 
In the present setting g £ BV(I); i.e., a function of bounded variation. It is well 
know that p £ := p £ m is a T e -acsm if and only if C £ p £ — p £ ; i.e., p e is a T e -invariant 
density. For each uj £ Q £ , we assume that T w satisfies the following conditions: 

(Bl) There exists a partition of 7, which consists of intervals {7j itJ }f =1 , 7^^(17^^ = 
for i ^ j, T itU := [a jU , c i+ i, w ] such that 

(i) for each i and for all < uj < e, e small enough, we have that (the quantity 6 
was introduced in the assumption (Al) above): 

[cifl + 8, Cj+1,0 - S] C [ci,u,Ci+i,u] C [ci t o - 5,c i+ i.Q + 5]; whenever c ifi ^ and 
Cj+1,0 7^ 1. In this way we have established a one-to-one correspondence between 
the unperturbed and the perturbed boundary points of Ii and If Ci.o and Ci iUJ 
are two such (left or right) corresponding points we will ask that Vi and Vw £ fl £ : 



(3.5) lim |cj o — Ci w | — (uniform collapsing of boundary points) 

e— >0 ' 

(ii) The map T u is locally injective over the closed intervals of class C 2 in 
their interiors, and expanding with inf^^ \T' u x\ > j3 > 2. Moreover, Vw £ Q, 
Vi = 1, • • • , g and Vx € [c;.o — <5, q + i.o + S] we have 

(3.6) lim \Tifi(x) — Ti tU (x)\ = (uniform convergence of maps) 
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(B2) Boundary condition: 

(i) if b <£ C , then T(b) = b and for all uj > 0, T u (b) = b; 

(ii) if b G C , then T(b-) < b < T{b+) and for all cu > 0, b G C u , where 

Cu) = { c i,uj}i=l- 

The last condition, as explained in section 2.4 of [23 does not allow the appear- 
ance of other holes, besides those around the infinitesimal holes, in the neighbor- 
hood of b. These conditions are satisfied when the noise is, for instance, additive: 
T u {x) — T(x) + uj. It is always well defined on the circle by taking the mod-1 un- 
folding; however, on the interval we should consider maps T for which the images 
|T w (a;)| lie in the unit interval for all u) G f2 e . In both cases the intervals of local 
monotonicity will be always the same and the differences between the perturbed 
and unperturbed images are uniformly bounded by s. Another example of a T u is 
shown in Figure 3. 



0.4- 



( Hr , 1 ^ 1 1 1 , 1 , } 

0.2 0.4 0.6 0.8 1 



Figure 3. A typical graph of one of the constituent maps T w of 
the metastable random map T e corresponding to the initial system 
shown in Figure 2. 

Under the above assumptions it is well known that (see for instance [5]) 

Proposition 3.2. For sufficiently small e > ; there exists a rj G (0, 1) and a 

B G (0, oo), such that for any f G BV{I) and n > 1, we have 

(3.7) ll^/IU<r? n ||/|| w + S||/||i, 

where • G {0,e}; i.e., both £ and C e satisfy a uniform Lasota-Yorke inequality as 
operators on BV(I). 



Among other things, Proposition 3.2 implies that 1 is an eigenvalue of C e and 
T E admits a stationary density p e . We further assume assume that: 

(Rl) p e is the unique invariant density of T e (ergodicity). 
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Remark 3.3. Since in this section we assumed tt £ to be finite, then condition (Rl) 
is satisfied if one of the maps T u has a unique acim. See Corollary 2 of [2"3"] . 
For instance, condition (Rl) is satisfied whenever the deterministic map T e is 
topologically mixing. 

3.3. Random holes in the perturbed system. We are interested in perturba- 
tions of T which produce "leakage" of mass from Ii to I r and vice versa. For this 
purpose, for each u G f2 £ we define the following sets: 

H Lu -iinT-^ir) 

and 

H r<u -.^irnT-^h). 

The sets Hi >UJ and H r>u , are called the "left hole" and the "right hole" , respectively, 
of the map T u . Thus, when T u allows leakage of mass from 1\ to I r , this leakage 
occurs when orbits of T u fall in the set Hi^- Similarly, when T w allows leakage of 
mass from I r to Ii , this leakage occurs when orbits of T w fall in the set H r 0J . 

For each of the left and right 'random open systems', we define a transfer operator 
which will be used to find the exponential escape rates from 7/ and I r . For this 
purpose we set 

and define for / £ L X (I;, 93(1;), f) 

(3.8) C L J(x):= f Ci fi (fl Xl J(x)de e , 

where Ci t o is the transfer operator associated with the map Tj J; := Ti : Ii — > J;. 

In a similar way, we can define a transfer operator C r ^ £ associated with the right 
'random open system'. 



3.4. Statement of the main result of section |3.3[ Using Proposition |3.2| and 
compactness arguments, it is well known that p e converges in the J^-norm to a 
convex combination of pi and p r . Our main goal in this part of the paper is to 
explicitly identify the weights in the convex combination. Since the acsm is ergodic, 
at least one of the two holes will have positive Lebsegue measure; we will suppose 
without loss of generality that the left hole will have this property. It will also carry 
positive mass for the left unperturbed invariant measure \xi since we assumed that 
the densities of fii and p r are positive in the neighborhood of the respective holes. 
We therefore define the limiting averaged holes ratio (l.a.h.r.) by 

(3.9) ^--'-HrlSI 

whenever the limit exists. We obtain the following theorem: 

Theorem 3.4. Let p £ be the unique % e -invariant density. Suppose the l.a.h.r exists 
then 

hm||p E -A)||i = 0, 
where po — api + (1 — a)p r and l.a.h.r = 
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Theorem 



2.1 



will give us additional information about the ratio jzr^ which ap- 



pears in Theorem |3.4| under the additional assumption: 



(B4) The map T is continuous on h* and Vcj, h* € H* >u . Moreover, (/*, /i*,T*), 
* € {I, r}, is mixing. 

Corollary 3.5. Assume condition (B4). Then 

(1) For sufficiently small e > 0, there exists an < e* i£ < 1, and a g* >e (x) > 0, 
with Jj g*^ £ dm — 1, such that 

(2) Moreover, 

,. l-e;, e a 
lim = . 

e->o 1 — e rj£ 1 — a 

Remark 3.6. As it was observed in (52], and this remains true in our case as well, 
the limit (3.9) surely exists whenever the perturbations open up holes whose size 
is first order in e. 

We also note that condition (B4) is satisfied, for instance, when the map T is 
continuous in the neighborhood of the infinitesimal holes and we perturb it with 
an additive noise. 

3.5. Technical lemmas and the proof of Theorem |3.4[ 

Lemma 3.7. We have J n (j, s (Hi >u )d9 s (u)) = j n n £ {H r ^)d9 £ {uj), where dji e 

f) til!) . 

Proof. We have 

(3.10) fiAIi) = [ tie{Ii)dO £ {u) = f fx e (X ltU) )de s (oj) + f tie(Hi,u)d9 £ (uj). 

Jn e Jn e Jn E 

Since T [ J 1 (/;) = Xi tU U H rjU , we also have 

(3.11) / /ieCC^MW = I He{Xi*,)d8c(u) + I f*e(H r , u )d6 e . 



Using (3.10), (3.11) and that [i £ is T e -invariant the lemma follows. □ 

Lemma 3.8. There exists a constant K , depending only on the map T , such that 
for < n < oo there exists a, depending on n, and e a , depending on a, such that 
Ve < e a and a < a, we have 

(1) £™1 is a C 1 -function on = [/i* — cr, /i* + cr], * G {I, r}; 

(2) for any x,y G Id, we have 

\L n e l{x)-L n e l{y)\<K-\x-y\; 

in particular, 

\C n e l(x)-C n e l{y)\<K.a. 

Proof. By condition (A4) we have for every n > 0, (TqCq) D Hq = 0. In particular, 
for a fixed n E N, we can choose cr' > so small so that, for < k < n, (Tq Co) D 
(h* — cr' , h* + a') = 0. We now prove that there is a possibly smaller subinterval 
(h* — cr,h* + cr), cr < a' upon which (1) and (2) of the lemma are satisfied. Let us 
give the explicit form of this operator. 
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The iterates of C £ are given by: 

= J ■ ■ ■ J C U1 (£ U2 (• • • £ w „ (1))) dfl e (wi) • • • d9 e {uj n ). 

Let us call . the ^-domain of injectivity of the map T LJj and call Tf~ u . the in- 
verse of T bJj restricted to Ii tU , . By Tip we denote the restriction of the unperturbed 
map to its Z-th interval of monotonicity; finally for simplicity of notation we drop 
the suffix Z, r from the intervals of the partition. We have: 

*«,„,... m (x) := (C Un _, (• • • C Ul (1))) (x) = 

\D(T Un o • • • o TuM^ o • • • o T^J( X ))\ x 1 ^::.'2>J x) - 

The sets 

(")ki,— ,k n . 

"ui,- ,U>„ ■ 

T kuui ° " ' ° T k n 1 . 1 ,u, n - 1 I k n ,^n n o • • ■ n T-^h,^ n 4 1>W1 

are intervals and they give a partition mod-0 of I = [0,1]; moreover the image 

,*„ . = T . . . T Qfcl,- ,fe„ 

= T Wn i kn<Un n T Un T UJn _ 1 i kuuJk _ 1 n • • • n T UJn T uln _ 1 o • • • o T Wl i kltU1 

for a given n-tuple {fc n , • ■ ■ , fci} is a connected interval. We will prove, by contradic- 
tion, that the function iUl (s) is smooth in a neighborhood of /i*, actually C 1 , 
when /i* is in the interior of one or several of the images H^l'V.'. ' t J^ described above. 
If V? Wn) ... , W i (#) is not smooth in a neighborhood of h* on any interval O D /i* , we can 
find a sequence (w n , • • • w n _;), 1 < I < n — 1 and at least an endpoint of one of the 
domain of injectivity of T Un _ l , call it c* u]n _ i , such that T UJn o ■ ■ ■ o T Un _, (c* ) € C 
We now proceed by induction and we begin by showing that is smooth in a 
neighborhood of /i* for any choice of u>. If not there will a point c; jW (see above, 
where we now mean that this point is one of the two boundaries of the interval 
Ii.u), such that T Ul ci^ £ {h* — a'/2, h* + c'/2). Let us now take the corresponding 
point cifl of T which will belongs to Co- We have 

\Ti,o(ci,o) - Ti u{ci,u)\ < \Ti,o(ci,o) - Tj i0 (cj w )| + |Tj i0 (cj w ) - T JiW (cj u)] 



By the uniform conditions (3.5) and ( |3.6[ ) we can take e smaller that a certain £ CT / 
depending on a 1 such that 

|7/,o(Q,o) - fj )0 (cj, w )| < sup |£>T;, ||ci,0 - Q,a,| < cr'/ 2 

and 

|T , ; ,o(q,c)-T ; , w ( C/!W )| <cr'/2 

and this implies 

|2j,o(cj,o) - < Cr' 

which contradicts the above statement T n Co n -f/o = 0. Fix a = a' /2; we now 
continue the induction process by supposing that for any sequence (uii, ■ ■ ■ the 
function $ w ... |U)1 is smooth on the interval (/i* — C, /i* + a) and we want to prove 
that the function ^ , tl) „ +lj ... is still smooth for the sequence (wj., • • • , w n +i). For 
that it will be enough to repeat the previous argument by noticing that 
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This shows that ^ Unj ... tUl (x) is smooth in a neighborhood of /i*. We now take 
two points x, y in the interval (h* — cr,h* + <r); we introduce the notation Tq := 



Tj, 1 , o • • • o T k u and we compute 



(3.12) 



E 



\D(T Un o • • • o T^)^ 1 ^))! \D{T Un o • • • o T^)^" 1 ^))! 



where we skip the characteristic function since both x, y will be in the same flu\\ 
We have 



D 



1 



(3.13) 



Since 



E 



■p-rn— 1 — fc rp 



yrn—l — k rp 



n 



3=0 



DT U , 



sup 



D 2 T u {x) 



< C\ < oo, 



DT u (x) 

which follows from our assumptions on the map T and its perturbations, and 

inf \DT u (x)\>P>2, 



the sum in (3.13) will be bounded by C\ times the sum of a geometric series of 
common ratio /3 : we call C2 the upper bound thus found. Therefore, by the 
mean value theorem, we have 



(3.121 <C 2 ]T \T^(x)-T^(y)\<\x-y\ ]T 



fei, 



\D(T ul c..oT Un )(T^(0)\ 



where C G (x, y). We can replace £ with a; (or y) by a standard distortion argument 
which works in our case by the assumptions we have on the maps T w . By modifying 
the constant Ci into a constant C3 which takes into account the distortion factor 
we finally have 



(3.14) 



(3.12) < C 3 |x- 



Vl k ^ kn \D(T Ul o...oT Un )(T^(x))\- 



We now integrate (3.14) over the noise and get, for any i,i/£ 



\C»l(x) - £™l(y)| < CsvCZl < <V(1 + B), 

where £™1 < 1 + B is obtained from the Lasota-Yorke inequality in Proposition 
3.2 The Lemma finally follows by choosing K = 63(1 + B). □ 



In [55], the authors dealt with the deterministic version of this section of our 
paper. In (25], the notion of the postcritical set of a map was used extensively to 
identify the location of discontinuities of the invariant density. However, in the 
random setting, the notion of a postcritical set of 1 e does not exist. Following the 
ideas of our Lemma |3.8| above, we develop an approach which is suitable for our 
random setting. We first recall a useful representation of a 1-dimensional function 
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of bounded variation^] For / £ BV(I), choose a version of / with regular discon- 
tinuities: for each x, f(x) = (lim y _ > . x - f(y) + lim J( _>. a ,+ f(y))/2. Then, / can be 
uniquely decomposed as / = f reg + f sal , where the regular term f reg is continuous 
and of bounded variation, with V(f re9 ) < V(f), while the singular (or saltus) part 
f sal is the sum of jumps 

r i = J2 s u h u , 

ues 

where S is a countable set, and 

{—1 if x < u 
— \ if x = u . 
if x > u 

This representation imposes the boundary condition f sal (l) = 0. Moreover, the 
variation of the singular part satisfies 



v(f sal ) = £ \ Su \ < v(f). 



Call 



n— 1 n—1 p 

n ' n ' Jo 

1 rt I n " lip 



fc=0 fc=0 ' 

where L% — £ Ul ■ ■ ■ C ulk . Note that, by (Rl), we have 

lim \\F n e - jO £ ||i = 0. 



We have already shown in Lemma 3.8 that the set of discontinuities of is given 

by 

where k n —j runs over the domains of injectivity of T Un _ . and dl denotes the end- 
point of the connected interval /. We rewrite this set as 

S n ,u> = U^TqSj^, where S jtQ := U kn _ j T UJn T UJn _ 1 ■ ■ ■ T Un _ j dI kn _ j ^ n _ j . 

We also write 

S n := U&S nt Q and Sj = U^o^a. 
Notice that Sj C S n . Finally, we write 

<S := U„>io„. 

Lemma 3.9. 

(1) S n is a finite set and S is a countable set. 

(2) [0, 1] \ S n is a finite collection of open intervals. 

(3) For any u> := (u) n , . . . , u>i), T a is C 2 on each open interval belonging to 
[0,1] \S n . 

Proof. Since f2 £ is a finite set, (1) follows by definition of S n and S. (2) is a 
consequence of (1). (3) follows by definition of <S„. □ 



The usefulness of such a representation in studying stability and response theory of 1- 
dimensional invariant densities was popularized by Baladi 5 . 
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Remark 3.10. As we have mentioned in footnote 4, if one relaxes the assump- 
tion that fl e is finite, even to the case where Q e is countable but not finite, then 
intractable complications arise. Indeed, if f2 e is countable but not finite, one can 
construct an example of a random metastable system such that in an interval, [u, v], 
say, each rational point is in the postcritical set of a different T u . For a system like 
this, the complement set [0, 1] \ S n becomes uncountable, and its intersection with 
[u, v] will not contain any intervaQ This makes it impossible to obtain regularity 
properties of p r e e9 , the regular part of the stationary density p £ . 

Lemma 3.11. The discontinuities of F n>£ belong to S n . 

Proof. We prove the lemma by contradiction. Namely, let a be a discontinuity 
point for F n s , and suppose that a (fc S n ,. Then Vw and < k < n— 1, we have that 



(3.15) 15m£*l(a:)=/:*l(a). 

Since, for sufficiently small e, each C u , uj £ f2 e , satisfies a Lasota-Yorke inequality 



with uniform constants (see Proposition 3.2) 



vc u f < nVf + b 

we concatenate the transfer operators and get the existence of a positive constant 
C, independent of the realization, for which 

H4i||=o<ll4ilUv<c. 



Therefore we can use (3.15) and apply the Lebesgue dominated convergence theo- 
rem to get 



71— 1 P 71— 1 



= lim - V / = lim F n , e (o;). 

x— >a n ' — ' /o a;— >a 

fe=0 

Thus, a contradiction. □ 



Let 



Lemma 3.12. 



#w = inf{j >l;u6 5,}. 



(1) The discontinuity set of p e is subset of S . If we write p s E al := J^ues s uH u , 
thenj: ueS ^ u>m \s u \<n-™C>. 

(2) p r e e9 is Lipschitz continuous. 



''One can easily construct a random system with such a bad behaviour. For example, let 
us write the formula of the first branch for the map To in Figure 2. We have To (a;) := 3x for 
< x < 1/6. For a fixed rational number e, Let fi e = Q n [0, e]. Then, each w £ r2 e is a rational 
number. Further, define the first branch of T ul {x) := (3 + 6u)x, for < x < 1/6. Then uj + 1/2 
is in the postcritical set of T u . In particular Si PI [1/2, 1/2 + e] = {uj + l/2}„gf2 e . Consequently, 
[0, 1] \ (S 1 ! n [1/2, 1/2 + e]) is an uncountable set and it does not contain any interval. 
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Proof. Since F n , E is a function of bounded variation, F n 



preg 



Fig. By Lemma 



3.11 



can be decomposed as 
the discontinuity set of F ni£ is contained in 
S n C S. Then we can write the saltus part of F n<£ as 

psal \ ^ 



ues 



iH», 



We first show that X)#«>m l s «,«l decays exponentially fast in m. 



To show this, 



observe that if m > n then Y" 1 _ e „ 



u>m \°u,n 



— and 



^ Kn| < V(F ntB ) < V(C k £ l) < C. 



Moreover, 



ueS; #«>0 



F„ 



1 



L F 



1 



l,e 



n n 

Now for #u > 1, i.e. u ^ Si, which means that Vw, each preimage T"T is continuous 
at u. Therefore the contributions to the jump s Uj „ come from the jumps of F n -i jS 
and are located at the points v such that T^v — u, precisely 

1 — / E \sv, n -x\dG e {J). 



n i] 



Then if < m < n we have 



E i 

#u>m 



< 



< 



< 



E - 

n-ll 



1 1 

V 



E i 



n rj 



E i s «,»-ii 

#u>m— 1 

^ ^ |^u,n— m| 
#tt>0 



< 



< 



C 



By using a diagonalization argument, we may find a subsequence rij such that for 
each u € 5, converges as rij — > oo to some number, which we call s' u . Let us 

construct the function 



psai 



ues 

It is clear that F^ £ converges in L 1 to F^ al and that for each m: J2 u es # M >m l s LI — 



r)~ m C. We now show that F^ ai coincides with p s e al and F reg , a limit point of i^ e /, 
coincides with p E es . Indeed, by (3) of Lemma 3.9 F ni£ is C 1 on the open intervals 
of [0, 1] \ S n . Therefore, on [0, 1] \ S n , the derivative of F ni£ coincides with the 
derivative of F^ eg (since F* a J is constant on each interval of [0,1] \ S n ). Then, 



using the arguments in Lemma 3.8 in particular estimates (3.131 and (3.14), we 



obtain a uniform estimate on the derivative of F^ eg on [0, 1] \S„. By (2) Lemma 
[0, 1] \ S n consists of a finite number of intervals. Thus, since by definition F^ eg 



3.9 



is a continuous function, we conclude that F^ eg is uniformly Lipschitz on [0,1]. 
Consequently, by the Arzela-Ascoli theorem, we can find a continuous function 
F reg such that some subsequence of {F^ 9 e } converges in L°° , and hence in L 1 , 
to F reg . By the uniqueness of the decomposition p e = p s s al + p 7 e e9 , we get that 

preg = p reg &nd psal = p sal Q 
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Lemma 3.13. 

(1) For any r > there exist an n, < n < CO, a a depending on n, and an 
e a , depending on a , such that Ve < e a and a < a, we have 

V IiP f < vT n C> < r, 

where I„ = [h* — &, h* + a], * € {I, r} and C" is the constant which appears 
in Lemma \3.1S\ 

(2) For any 1,1/6 Ig., there exists a uniform constant K such that 

\p^(x)-p^(y)\<K-\x-y\. 



Proof. Choose n large enough so that rj~ n C < r. Then, using (1) of Lemma 3.8 for 
n > such that r]~ n C < r we can find a a, depending on n, and an s a , depending 
on a, such that Ve < e a and a < a, £™1 is a C 1 -function on 1„. Consequently, 



using Lemma 3.12 with m = n, we obtain part (1) of the lemma. Part (2) of the 
lemma is a consequence of part (2) of Lemma 3.12 □ 



We now observe that by the Lasota-Yorke inequality of Proposition 3.2 and the 
compactness of the unit closed ball of || • \\bv in L\ there are subsequences of 
values e£ !1 E for which converges in the L 1 norm to the density of an acim on 
/. Obviously, such a density will surely be a convex combination of pi and p r . We 
would like to show that those accumulation points are always the same so that p £ 
will admit a limit when e — > 0. The first step is to characterize the coefficients of the 
convex combination in terms of the behavior of the density in the neighborhoods 
of the holes. We first have: 

Lemma 3.14. Let po be the accumulation point of the convergent subsequence pg 
with i € Q e ; then there exists < a < 1 such that^ 

(1) po = apt + (1 - a)p r ; 

(2) lira sup sup \pg(x) - api{x)\ = 0; 

(3) lim sup sup \pg(x) — (1 — a)p r {x)\ = 0. 

Proof. We argued above that (1) is true. The proof of (3) will be identical to that 
of (2). Thus, we only prove (2). Observe that for all u £ !! E - we have Hi^ C H^g. 
Therefore, 

sup sup \pg(x) — api(x)\ < sup \pg{x) — api(x)\. 
Thus, it is enough to prove that 



(3.16) 



lim sup \pg(x) - api(x)\ = 0. 



We will prove (3.16) by contradiction. For this purpose we will suppose that the 
there exists a C > and a subsequence e' — > of e such that for each e' there is 
an x e i £ FIi,£', x e > — > hi with \p e i{x e i) — api(x e >)\ > C. By Lemma 1 of [55], for 
sufficiently small a > 0, for any x € [hi — a, hi + a] := I s and x s i € I s we have 



(3.17) 



\api(x e >) -api(x)\ < 2C/5. 



^Later in the proof of (1) of Theorem 
show that l.a.h.r = , - . 

l — a 



3,1 



we will find 



explicitly. In particular, we will 



W. Bahsoun And S. Vaicnti 



2:S 



Choose n big enough so that C'rj ™ < C/5. By Lemma 3.13 we can find a a 
depending on rt, and an e^, depending on er, such that Ve' < e' a and a < a, we have 

Vj,p s f < C/5, 

and for any x, y € ig-, there exists a uniform constant if such that 

\p™°{x)-pl«°{y)\<K-\x-y\<Kd. 

We make sure to take a small enough such that Ka < C/5, and a < er. Then for 
x, x £ > £ 1 5 we have 

(3.18) M^v) - ^(z)] < - p7\*)\ + \p7 9 (^) - p r E f 9 (x)\ < 2C/5 



Thus, for any x € Is and x e i <E Ig., by (3.17), (3.18) and the assumption that 
\p e >{x s >) - api(x e ')\ > C, we obtain 

(3.19) 

\p e '{x) - api(x)\ > \p £ >(xs>) - api(x £ ,)\ - \api(x) - api{x e >)\ - \p e '(x £ ') - p E '(x)\ 
> C-2C/5-2C/5 = C/5. 



(3.16 1 holds and t 



The estimate in (3.19) contradicts the fact that lim||pj — P0II1 = 0. Therefore, 
ic lemma follows. 



□ 



We are now ready to prove Theorem |3. 4 



Proof, (of Theorem 3.4) Let po = api + (1 — a)p r be the limit of p £ obtained in 
(1) of Lemma 3.14 We will show that a/(l — a) = l.a.h.r. This will imply the first 



part of Theorem |3.4| We first observe that 
(3.20) 



p,e(H ltU )d9 g (u}) 



O f J H. 



m s j H-, 

apidxd9i(u>) - 



pgdxd9 £ (uj) 

(pi - api)dxd9i(u}) 



n g Jh 



= a pi{H l ^)d9 i {uj) + O sup sup \p £ (x) - api(x)\ \ I m(H liUJ )d9 e -(u). 

Let hi be the infinitesimal hole in //. By condition (A4) po is continuous on hi. 
This implies lim^o mfg;' " ) = Plfoi) > by condition (A5). Therefore, by using 
(2) Lemma 3.14 and (3.20), we obtain 



(3.21) / fie{Hi^)d6e(u) = a / ^(H t ^)d6 e (u) + o(l) • / pi{H hu )d9 s {oj). 
Jsn i Jn € Jn € 

Similarly we can obtain 

(3.22) f iig(H r , w )d$e(w) = (1 - a) [ p r {H r ,„)dO e {u) + o(l) • f p, r (H r ,„)d0g{u). 

JCli: JQi: JCli: 



Using Lemma 3.7 together with equations (3.21) and (3.22) complete the proof of 
Theorem KM □ 
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4. DISCUSSION 

In this paper we studied both open and metastable systems. In Theorem |2.1| 
we rigorously derived first-order asymptotic escape rate formulae for piecewise ex- 
panding maps of the interval with random holes. In |27| the authors studied a 
deterministic sequence of nested holes and derived a first-order asymptotic escape 
rate formulae for piecewise expanding maps of the interval with holes. In this pa- 
per the sequence of holes is random. Moreover, the random holes we consider are 
not necessarily nested. Consequently our work generalizes the pioneering work of 
[27] to the random setting. It would be interesting to explore whether one can 
obtain higher order approximations in the random setting, even for a simpler class 
of maps, such as circle maps. A study for higher order approximations of escape 
rates for deterministic holes using circle maps was carried in [2]. Other possible 
generalizations, such as higher dimensional systems, were discussed earlier in the 
introduction. 

In the second part of the paper we studied the problem of approximating the 
unique stationary density for a random system which initially contains two separate 
ergodic components. We have shown in Theorem |3.4| that the unique stationary 
density of the random system can be approximated by a weighted combination of 
the two initially separate ergodic components of the deterministic system. Further- 
more, we have shown that the weights are given by the ratio of the escape rates 
on the individual ergodic components. A deterministic perturbation version of the 
above result was first studied in [221 for piecewise expanding maps. A version of 
|22j for intermittent maps was obtained in [3J. It is worth restating (see introduc- 
tion) that it would be also interesting to obtain a random version in the case of 
intermittent maps. 

Acknowledgment We would like to thank anonymous referees for their sugges- 
tions which improved the presentation of the paper. 
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